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ABSTRACT 

Using the CLEO-II detector at the CESR e+e“ storage ring running at the 
center-of-mass energy around 10.6 GeV, we have studied exclusive production 
of charmonium states in two-photon collisions. Employing a dataset comprising 
3.0 fb“^ we have searched for decays of the rjc, XcO and Xc 2 mesons resulting in 
four charged hadrons. We report on the measurements of the cross sections and 
two-photon widths of these charmonium states. 

1. Introduction 

In this paper we present the results of a study0 of charmonium states exclu¬ 
sively produced in two-photon collisions. The processes considered here are e’''e“ ^ 
e'^e~Rcc, where Rcc is either rjc, Xco or Xc 2 mesons. We measure their production 
rates when neither the scattered electron nor the positron is detected (“untagged” 
mode). In these processes a meson is coupled to two space-like photons one emitted 
by the electron, the other by the positron. From the measured production rates we 
obtain two-photon widths of the mesons. The following decay channels are used in 
this study: rjc — >■ Xco —*■ 7r+7r“7r+7r“ and Xc2 —*■ 7r+7r“7r+7r“. We also search 

for the rjc decaying into four charged pions. 


2. The Detector 


A brief description of the major components of the CLEO-II detectoiH and its trig¬ 
ger system! can be found in these proceedings.0 Additional features of the detector 
relevant to this analysis are: Charged particle identification provided by measure¬ 
ments of ionization losses in the drift chamber (dE/dx) and time-of-flight measure¬ 
ments from scintillation counters; Electron identihcation making use of the dE/dx 
and E/p ratio, where E is the energy of the cluster produced by the electron can- 
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didate in the calorimeter and p is the magnitude of the charged track momentum 
measured using the drift chambers; Muon identihcation uses the depth of particle 
penetration into the hadron hlter. This analysis utilizes data recorded with triggers 
which are designed to be efficient for events containing multiple charged tracks. 

3. The Event Selection 

In order to select event candidates we require that the following selection criteria 
be satisfied. 

• There are exactly four charged tracks with zero net charge. 

• The root-mean-squared residual of the hits associated with each track is less 
than 0.4 mm. 

• The total detected energy is less than 6.0 GeV. 

• The total amount of energy collected in extra cluster(s)@ not associated with 

any of the charge tracks is less than 500 MeV for the mode and less 

than 650 MeV for the four charged pions mode. 

• There are no particles identified as electrons, positrons or muons. 

Additional criteria specific to the considered decay channel are discussed in the next 
two sections. 

4. The ?7c —> KgK^TT^ Decay Channel 

The decay —>• tt+tt” is identified taking advantage of the large decay 

length (cr = 2.675 cm). We look for pairs of oppositely charged tracks intersecting 
at a transverse distance (he. projected distance in the plane perpendicular to the 
beam axis) greater than 1 mm from the primary interaction point. The momentum 
of each pair is calculated using the track momenta evaluated at the secondary vertex. 
The energy is calculated assuming that tracks are produced by pions. The invariant 
mass of each pair is then calculated. The resulting mass spectrum in Fig.^ shows 
clear evidence of decays I A pair is accepted as a candidate if its invariant 
mass falls within 12.0 MeV/c^ of the peak position. We reject events if the two 
remaining tracks are consistent with a decay: he. if their invariant mass falls 
within 12.0 MeV/c^ of the peak position of the previously discussed distribution. 
To supress beam-gas background events and to select events of reasonable quality we 
require that the two tracks recoiling against the come from the nominal interaction 

“The position of the peak in the mass distribution is obtained by fitting this distribution with a 
Breit-Wigner function for the signal and a constant for the background. 
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region: their closest approach to the primary interaction point is less than 50 mm 
in the longitudinal direction and less than 5 mm in the transverse plane. These 
tracks are assigned the K and tt masses and the mass value is calculated 

for each of the two possible mass assignments. The upper tail probability of the 



Fig. 1. a) The invariant mass spectrum for the candidates in four-prong events. The curve shows 
the result of the fit described in the text, b) The transverse momentum (p_l) distribution for the 
candidates with the invariant mass from 2.7 GeVjt? to 3.3 GeVjc^. 


distributions (Putp) is calculated for each mass assignment using TOF and dE/dx 
information. In the ideal case, the correct choice of particle combinations produces 
a uniform P^^tp distribution while incorrect particle combinations tend to congregate 
near zero. To supress unwanted background events which do not contain charged 
kaons Putp is required to be greater than 0.05. Most candidates fail this condition 
while the remaining events are dominated by a single acceptable mass assignment. 
Events are given a net weight of one. For events with two acceptable combinations, 
each possibility is given a weight proportional to its value of Putp- 

In the two-photon untagged events substantial energy and momentum are carried 
away by the electron and positron scattered at small angles. However, since the 
transverse momenta of scattered electrons are small, that of the hadronic system 
should also be small. Hence, the distribution of the transverse momentum of the 
hadronic system {p±) peaks at zero as shown in Fig.|^. Events with a transverse 
momentum {p±) of the hadronic system less than 200 MeV/c are accepted into the 
hnal invariant mass plot for the candidates which is shown in Fig.|^a. To 

obtain the number of events in the signal peak we £t this distribution with a signal 
shape function 0 obtained by running a GEANT-based detector simulation program. 
A smooth background contribution is approximated by a power function of the form 
N{W) = NqW^, where Nq is a normalization parameter, x is a power parameter 


“The signal shape function is a combination of two Gaussians with different widths and the same 
peak position. It reflects the variation of the detector resolution over the angular and energy ranges 
of the decay products. 
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and W is the invariant mass. The nnmber of observed events is shown in Table |^. 
5. The TT+TT TT+TT Decay Channel 


To identify the decay mode into fonr charged pions we reqnire that all charged 
tracks originate from the interaction region. To satisfy this reqnirement each track 



Fig. 2. The mass distributions for the candidates: a) b) tt+tt tt+tt for tJc mass range 

and c) 7r+7r“7r+7r“ for Xc mass range. The curves show the results of the fits described in the text. 


must have a separation from the beam crossing in the longitudinal direction of less 
than 50 mm and in the transverse plane of less than 5 mm. The upper tail probability 
of the distribution is determined from TOF and dE/dx measurements assuming 
all four tracks to be pions. Candidates with Putp greater than 0.1 are accepted into 
the 7r+7r“7r+7r“ invariant mass plot. In addition, we require that px of the hadronic 
system be less than 200 MeV/c. Signals for both the Xco and Xc 2 are seen in Fig.|^c. 
The numbers of signal events are obtained in the same way as in the previous section. 
The background contribution is approximated by an exponential function. 

We do not observe rjc in the decay mode into four charged pions. The distribution 
of the invariant mass around the region where one expects pc is shown in Fig.^j. The 
vertical arrow shows where an pc signal is expected. 

6. Results 

We obtain the detection efficiency by running a detector simulation program which 
includes trigger simulation. The trigger is efficient for 70% of all events and for 94% of 
events which would satisfy analyses selection criteria. These efficiencies agree within 
5% with trigger efficiencies measured in data using redundancies between different 
triggers. For simulation of the two-photon processes we use Monte Carlo generators 
based on the formalism of Budnev.i The intrinsic properties of the charmonium 
resonances are described by a double pole form factor model.0 The pole mass has 
been varied from zero to the J/'ip mass. We estimate that the uncertainty in the pole 
mass and form factor model are sources of a relative systematic error of 10%. Other 
sources of systematic errors include uncertainties in the detection efficiency due to 
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Table 1. Summary of rjc, Xco and Xc 2 measurements. 


Channel 

Number 

observed 

Branching 
Fractions (%)i 

Br • F..y^ 

(0.01 X KeV) 

r 

(KeV) 


54.1 ± 12.6 

1.5 ±0.4 

6.5 ± 1.5 ± 1.1 

4.3 ±1.0 ±0.7 ±1.4 

i 

+ 

1 


1.2 ±0.4 

< 2.7 ( 90% CL ) 

< 2.4 ( 90% CL ) 

+ 

1 

47.2 ± 15.5 

3.7 ±0.7 

6.4 ±2.1 ±1.5 

1.7 ±0.6 ±0.4 ±0.3 

Xc2 —^ 
7r’'“7r“7r’'“7r“ 

41.9 ± 13.8 

2.2 ±0.5 

1.5 ±0.5 ±0.3 

0.7 ±0.2 ±0.1 ±0.2 


the event selection criteria, trigger and detector simulation and tracking. Overall 
detection efficiencies are 12 %, 20 % and 21 % for the rjc Xco 47 r^ and 

Xc 2 decay channels, respectively. 

The results of this study are summarized in Table |l|. The first two errors quoted 
are statistical and systematic. The third error in the two-photon widths arises from 
the uncertainties in the corresponding branching fractions. Our measurement of the 
T^^{ric) = 4.3 KeV using the rjc decay mode is consistent with previous 

measurements which are in the 5-28 KeV range, though most recent and more precise 
measurements fall below 10 KeV. The upper limit on r.y.y (r^c) derived from the decay 
rjc 7 r+ 7 r“ 7 r+ 7 r“ is T^^{ric) < 2.4 KeV at the 90% confidence level. Our measurements 
of r.y.y for Xco and Xc 2 mesons are consistent with results of previous experiments.! 

7. Discussion and Summary 

The nonrelativistic quark model predicts that the two-photon width of the bound 
qq state with angular momentum I is proportional to the square of the l-th derivative 
of the qq pair wavefunction at the origin {i.e. at contact). To calculate the two-photon 
width one needs to know this wave function, or in terms of the Fourier transform, the 
momentum distribution of the qq pair inside the meson. This momentum distribution 
can be obtained from the nonrelativistic potential model. However, due to this non- 
perturbative calculation the overall normalization of the width is not very reliable. 
This uncertainty can be partially removed when one considers the ratio of two partial 
widths of the same meson or similar mesons. The ratio of the widths may be calculated 
more accurately using a relativistic approach. 

For the rjc meson a nonrelativistic model! with next-to-leading order quantum 
chromodynamics corrections predicts: /T = (4/3)(l -|- l.Qbas/vr). 

Using the measured width T{J/ip ~ 5.3 KeVS with a choice of the Ug ~ 

0.28, this yields T{ric 77 ) ~ 8.2 KeV. This calculation assumes that the rjc and 
J/ip wavefunctions are identical at the origin, ignoring hyperfine mass splitting and 
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coupled-channel effects. Indeed, the hyperfine effects would be likely to enhance 
the wavefunction of the relative to that of the Some evidence that the 

wavefunctions are not identical comes from the supression of the J/xjj ^ 'yrjc magnetic 
transition. On the other hand a description which incorporates relativistic effectJl^ 
with a wavefunction obtained from a potential model predicts T^^{r]c) ~ 4.8 KeV. ^ 
This calculation agrees with our measurement of the r^^{T]c) = 4.3 ± 1.0 ± 0.7 ± 1.4 
KeV obtained in the mode. 

For the Xco and Xc 2 states a naive nonrelativistic quark model predicts r.y.^(xco) ^ 
T^^{Xc 2 ) = 15/4 : 1. The next-to-leading order quantum chromodynamics corrections! 
change this ratio by a factor of (1 -|- 0.2as/7r)/(l — lOa^/(Svr)). By evaluating this ex¬ 
pression at as ~ 0.28 we obtain T^^{xco)'y-y{Xc 2 ) ~ 7.3. The relativistic calculation^ 
predicts this number to be around 2.8. Only the latter value is consistent with our 
measurement of 2.4 ± 1.1 ± 0.7. 
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“This is evaluated assuming the mass of the c quark (wc) to be 1.4 GeV. A somewhat larger rric 
= 1.6 GeV leads to the r.y.y(? 7 c) = 3.4 KeV. The choice of the c quark mass becomes even more 
important for P-wave qq states because of the higher power of the c quark mass involved. 
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